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Abstract. The highest-mass stars have the lowest frequency in the stellar IMF, and 
they are also the most easily observed stars. Thus, the counting statistics for OB stars 
provide important tests for the fundamental nature and quantitative parameters of the 
IMF. We first examine some local statistics for the stellar upper-mass limit itself. Then, 
we examine the parameter space and statistics for extremely sparse clusters that contain 
OB stars, in the SMC. We find that thus far, these locally observed counting statistics are 
consistent with a constant stellar upper-mass limit. The sparse OB star clusters easily 
fall within the parameter space of Monte Carlo simulations of cluster populations. If the 
observed objects are representative of their cluster birth masses, their existence implies 
that the maximum stellar mass is largely independent of the parent cluster mass. 



1. Introduction 

Since massive stars are exceedingly rare, any local, empirical determination of their 
statistical properties is necessarily a stochastic problem. This applies to determining 
the slope of the initial mass function (IMF) and stellar upper-mass limit. While the 
properties of massive stars are inferred from models for both stellar atmospheres and 
stellar evolution, we can begin to set constraints on the statistical properties of the pop- 
ulation once we believe we can make consistent, systematic estimates of the stellar 
properties. However, it is important to bear in mind the systematic uncertainties inher- 
ent in these determinations (Massey, these Proceedings). Here, I will examine a few 
stochastic analyses of local observations of the OB star population. 



2. UP: The Stellar Upper-Mass Limit 

If the IMF truly behaves as a probability density function, which is indeed the way that 
we ordinarily assume that it does, then in principle we can collect data from many differ- 
ent clusters and combine them to increase the size of our sample to beat d own stochastic 



noise . We applied this technique to evaluate the stellar upper-mass limit (l Oey & Clarke 



2005h . using data from local OB associations. iMassey et all (1 19951) obtained uniform. 



high-quality spectroscopic classifications of the upper IMF in a number of OB associ- 
ations in the Milky Way and Magellanic Clouds. To evaluate the upper-mass limit, we 
are interested only in clusters having ages < 3 Myr, which are young enough so that we 
are observing these systems before the most massive stars have expired as supernovae 
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Figure 1. Expected maximum mass (m,nax) vs upper- mass limit niup of the parent 
distribution, for ensembles of N, = 100, 250, and 1000 stars having masses > lOM©. 
Dotted lines assume the parent IMF has the form of a truncated Salpeter power law, 
and solid lines assume the form in equation |2] The dashed line shows the locus for 

<'«max) = '«up- 



(SNe). Eight objects studied by Massey et al. (1995) qualify: IC 1805, Berkeley 86, 
NGC 7380, NGC 1893, NGC 2244, Tr 14/16, LH 10, and LH 117/118. These clusters 
contain, respectively, 24, 10, 11, 19, 12, 82, 65, and 40 stars having masses > 10 Mq, 
for a total of 263 stars. In addition to this sample, Massey & Hunter (1998) used the 
same techniques to obtain the census of the upper IMF in the R136a super star cluster 
in the 30 Doradus complex in the LMC, using HST/STIS spectroscopy. Their strong 
lower limit for the number of stars having masses > 10 M© is 650. Combining the 
samples of ordinary OB associations and R136a yields a grand total of 913 stars. 

The maximum stellar mass in the entire sample is around 120 - 150 M^. We can 
compare this with the average expected maximum mass for an IMF with a 'Salpeteij 
([1955) power-law slope y = 2.35, for an assumed upper-mass limit Map of the parent 
distribution, and a given number of stars A^* in the ensemble. As before, A^* includes 
only stars of mass m > 10 Mq. The expected maximum mass (nimax) is given by: 



(mmax) = m 



up 



Jo iJo 



(l>{m) dm 



N 



dM 



(1) 



Figure [T] shows the average expected (mmax) as a function of mup for three different 
values of A^* (dotted lines). For a true upper-mass limit m^p = 1000 Mq, an ensemble 
having A^* = 250 has an expected (mmax) ~ 450 Mq; an ensemble having A^* = 1000 has 
("^max) ~ 650 Mq. These values of A^* are similar to those for our observed combined 
samples of Milky way and LMC associations, excluding and including R136a. We note 
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that the values of (mmax) do have some dependence on the assumed form of the parent 
IMF. The dotted lines in Figure[T]show the results from equation [T]for a simple Salpeter 
power-law truncated at Wup. However, if we adopt a "softer" truncation of the form. 



n{m) dm oc 
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then instead equation [T] yields the solid lines in Figure [TJ 

In any case, as far as we know, no stars approaching these expected m^ax values 
on the order of a few hundred M© are known to be observed, and the example calcu- 
lations assume a true m^^ of merely 1000 M©, much less infinity. Indeed, inverting 
the argument, the observed valu es of mmax ~ 150 M© imply that the parent niup has a 
similar value. Elmegreen' (2000') used this reasoning to estimate that if m^^ = oo, then 
somewhere in the Milky Way, there should be a star having m^ax = 10,000 M©. No star 
remotely approaching this mass is suggested to have been seen. 

However, given that the preceding is based on only 9 clusters, we may worry that, 
because of stochasticity, these may not be representative of the massive star popula- 
tion. We can quantify the degree to which we might be this unlucky by evaluating 
the probabilities pim^ax) of obtaining the observed m^ax seen in each cluster, for an 
assumed wiup. Under normal circumstances, we should be uniformly lucky and un- 
lucky, and so these pim^ax) should represent a uniform distribution when mup corre- 
sponds to its actual value in the parent distribution. Figure |2] shows the histograms 
of /j(mniax) for assumed wiup = lO'*, 200, 150, and 120 M©. The probabilities P that 
these histograms are drawn from uniform distributions for these respective cases are 
P < 0.002, < 0.02, 0.12, and < 0.47. Thus we confirm that, only when m^p has values 
similar to the observed mmax> do we see significance in P, demonstrating that mup in this 
part icular samp l e is in deed around 150 M© at the significance implied by P. We note 
that lAban et all (l2006l) point out that the statistical estimator for the maximum value of 
a truncated, inverse power-law distribution is indeed the maximum value in the dataset, 
corresponding to the observed m^ax- 

It is especially remarkable that this apparent value of ~ 150 M© is seen over 
a great range in star-forming conditions. It holds for both ordinary OB associations in 
both the Milky Way, a massive, spiral galaxy; and the LMC, a lower-mass, irregular 
galaxy with presumably a lower-pressure ISM. At the same time, the same upper-mass 
Umit applies in the R136a super star cluster, an extreme environment in the LMC. And 
IFiger (2005) finds the same upper-mass limit in the Arches cluster near the Galactic 
Center, which is another extreme, yet very different, environment. We caution that 
these results are based on the inferred stellar masses for these systems, taken at face 
valu e, and assuming that the highest-mass stars have not yet expired. 

Koen (120061 ) used the same dataset for R136a (iMassey & Hunted 1 19981) to evalu- 



ate mup and the IMF slope y from the cumulative stellar mass function. The cumulative 
distribution more clearly reveals the existence of a truncated mup. Koen simultaneously 
fit mup and y for both sets of masses provided by Massey & Hunter, based on two dif- 
ferent calibrations for spectral type to stellar effective temperature. The fit was carried 
out with both a least-squares method and a maximum likelihood method. Across all 
four cases, the results again are largely consistent with a Salpeter slope and mup ~ 150 
Mq. Fixing mup - oo results in a very steep slope, y ~ 4 to 5, strongly inconsistent 
with the Salpeter value. An infinite mup is eliminated with a probability of < 0.0005 
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Figure 2. Histograms of pim^^^) for assumed - 10'*, 200, 150 and 120 M©, 
as shown. For unbiased conditions, piniitax) should show a uniform distribution. 



and < 0.002 for the least-squares and maximum-likelihood methods, respectively, con- 
sistent with our results (Oey & Clarke 2005). 



3. Sparse OB Star Groups 

Given the apparent robustness of the upper-mass limit in a variety of local environ- 
ments, we seek to examine other extreme circumstances. An especially interesting 
case is field OB stars. While a substantial fraction of these are likely to be runaway 
stars ejected from clusters, many may well be members of small, low-mass clusters. 
Oey et alj (120041) demonstrated that individual SMC field OB stars, as defined by a 
friends-of-friends algorithm, do fall smoothly on a continuous power-law distribution 
in the number of stars per cluster A^*, which is akin to the cluster mass function. This 
supports the scenario that most field OB stars are members of low-mass clusters. 

We obtained F555W and F814W SNAP observations of eight, appare ntly isolated 



OB s tars from among these SMC field stars, using the HST ACS camera (iLamb et al. 



20101) . Applying both a stellar density analysis and a friends-of-friends algorithm, we 
confirmed that three of these objects have stellar density enhancements, corresponding 
to sparse stellar groups. An additional object registered a positive signal for compan- 
ions using the friends-of-friends algorithm, but not the stellar density analysis. The 
remaining four stars appear isolated to within the detection limit of F814W = 22 mag. 
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Radial velocity observations subsequently revealed two of these stars to be runaway 
stars. This leaves two stars remaining as candidates for in-situ isolated, field OB stars. 

We constructed Monte-Carlo simulations of cluster populations to examine the 
parameter space occupied by our observed, sparse OB star groups. We first adopted a 
cluster mass function (MF) given by a simple power-law distribution: 

A^(Mci) dM^i oc M-^ dMd , (3) 

where N{Mci) is the number of clusters in the mass range to + dM^. Each 
cluster drawn from this distribution is then populated with a stellar IMF drawn from a 



Kroupal (120011) IMF, which has the form. 



, , J f dm , O.O8M0 < m < 0.5 Mq 

nun) dm < 7 ^5 , ^1., ^^r^i., (4) 

^ ^ [ m'^-^^ dm , O.SMq < m < 150 Mq ^ ' 

We adopt a default model having a power-law index of -2 for the cluster MF (equa- 
tion O, and lower-mass limit of Mdjo - 20 M©. This model best reproduces the 
frequency of single-0 star clusters in the SMC and Milky Way (see Lamb et al. 2010). 

Figure [3]2 shows the mass ratio distribution of the second-highest to highest mass 
stars mniax,2/"^max vs m^ax in these simulations. Only clusters having a single OB 
star, defined in this figure as having mass m > 18 Mq are shown. Our observed ob- 
jects are overplotted with the black squares, showing upper limits for the apparently 
isolated objects. We see that while our observed objects appear to fall in a densely- 
populated region of the parameter space, the simulations peter out at the lowest values 
of 

'Wmax,2/'Wmax» sincc the probability of drawing an extremely low mass ratio is tiny. 
All of our data fall within the lowest 20th percentile in ffimax.a/'Wmax. and these frequen- 
cies also result when our clusters are simulated based on cluster membership number 
Nt instead of a MF (Figure[3^). This model analogously uses the form A^7^> the Kroupa 
IMF, and a lower limit of A^*jo - 40 stars, which is equivalent to Mdjo used above. 
That our objects all fall in this low-frequency regime is partly due to our selection of 
apparently isolated OB stars as the HST targets. We note that the simulations reproduce 
this regime fairly easily. 

Figure |4] shows the relation between the m^ax and M^x for the simulations based 
on the cluster mass function (equation O, now showing all clusters having at least one 
OB star. Since there is no relation imposed between the mmax and M^, the simulations 
occupy a large parameter space. The solid lines show the 10th, 25th, 50th, 75th, and 
90th percentiles of m^ax for a given Md, and the dashed line shows the mean, equiv- 
alent to equation [J Our observed objects are shown with the black squares as before, 
with Mci computed by assuming that a fully populated IMF exists below the detection 
thr eshold of ~ 1.5 Mp. Open diamonds show the sample of observed clusters compiled 
bv lWeidner eTaP (1201 Ot) . 

We see that for almost all of our observed objects, 90% of the simulated clusters 
have lower nimax for a given Md, whereas the sample from Weidner et al. (2010) is 
largely at the opposite extreme. Our data, taken at face value, are in substantial conflict 
with the premise of a well-defined m^^x - ^ci relation, proposed by Weidner & Kroupal 
(12005.) . The observed data show a scatter spanning two orders of magnitude in Md for 
the wTmax in our sparse OB groups. While it is essential to understand the consequences 
of cluster dynamical evolution, this range in values, which is also consistent with the 
unconstrained simulations, strongly suggests that a mmax - ^ci relation is weak at best. 
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Figure 3. Mass ratio of the two highest-mass stars vs m^ax for cluster simulations 
based on a cluster MF (panel a; top) and on a cluster membership function (panel 
b; bottom). Only simulated clusters having a single OB star are shown, with gray 
points. Black squares correspond to our observed sparse O star groups, with three 
objects showing upper limits. 
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Figure 4. ?Wn,ax vs Md in simulations based on a cluster MF, for all objects 
with one or m ore OB stars. Our observed objects and the sample compiled by 
IWeidner etaP 12010) are shown in black squares and open diamonds, respectively. 
Solid lines show the 10th, 25th, 50th, 75th, and 90th percentiles of mmax, and the 
dashed line shows the mean. 



as also found by iMaschberger & Clarkd (l2008h . iTesti et al. 1 (119991) discovered sparse 
groups around Galactic Herbig Ae/Be stars, which are newly formed objects, and thus 
their parent systems are unlikely to be strongly depleted by dynamical evaporation. Our 
observations extend such findings to OB stars. 

The existence of a relation between m^ax and is widely debated. A clear re- 
lation would strongly impact the integrated galaxy IMF (IGIMF; Weidner & Kroupa 
2005), affecting interpretations of stellar populations, star-formation histories, and in- 
ferred galaxy evolution. Moreover, competitive accretion theories for star formation 
predict the existence of such a relation: wimax °^ 

mII^ (iBonnell et al.ll2004l) . It is there- 
fore of great interest to further investigate whether sparse OB star groups are represen- 
tative of their birth conditions. 



4. Summary 

With samples of uniformly derived OB star masses, we can begin to quantitatively eval- 
uate the properties of the stellar upper IMF, based on a stochastic understanding of its 
nature. We demonstrated the existence of a stellar upper-mass limit niup ~ 150 M© for 
a Salpeter slope, using data for the massive star census in a sample of ordinary Milky 
Way and LMC OB associations, and the super star cluster R136a in 30 Doradus (Oey 
& Clarke 2005). Based the probability of observing the highest-mass stars in their re- 
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spective clusters, we confirmed the existence of a strong deficit above this value of wiup. 
Koen (2006) examined the cumulative distribution function of the stellar masses, and 
by jointly fitting m^^p together with a power-law distribution for R136a, he confirmed 
a value for m^p ~ 150 and a Salpeter-like slope as the best-fit values for these pa- 
rameter estimations. The variety of environments in which this value of the upper-mass 
limit applies is remarkable: ordinary OB associations, the R136a super star cluster, and 
the Galactic Center environment. 

To further examine the robustness of this upper-mass limit, we searched for low- 
mass companions around eight SMC field OB stars with HST/ACS. We confirmed the 
existence of sparse groups associated with 3 - 4 of these field massive stars, while 2 
targets are runaway OB stars, and 2-3 remain candidates for isolated OB stars that 
formed in situ (Lamb et al. 2010). We generate Monte Carlo simulations of cluster 
populations based on simple power-law sampling for both cluster and stellar masses, 
and we find that the cluster lower-mass limit is most consistent with the frequency of 
SMC and Galactic field O stars at Mdjo ~ 20 M© or A^*jo ~ 40 for a Kroupa IMF. 
Our sparse OB groups all fall in the lowest 20th percentile in minax,2Mmax5 regardless 
of whether the clusters are populated by M^i or by A^*. 

Our sparse OB groups generally fall in the highest 10th percentile of mmax for a 
given Mci, at an opposite extreme from the sample of objects compiled by Weidner et 
al. (2010), which fall in the locus of the most massive objects, at a given mmax, in our 
simulations. Taken at face value, our results therefore contradict the existence of a well- 
defined relation between wzmax and Mc\. If our sample is not predominantly the result of 
dynamical evaporation, then this finding may pose difficulties for the predicted steep- 
ening of the IGIMF based on a suggested Wmax - Mc\ relation (e.g., Weidner & Kroupa 
2005) and competitive accretion theories for star formation (Bonnell et al. 2004), which 
also rely on the existence of such a relation. 
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